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CONTINUOUS ON-LINE TITRATIONS 
BY FEEDBACK BASED FLOW RATIOMETRY 

BACKGROUND OF THE INVENTION 

[001] This application claims the benefit of U.S. Provisional Application No. 

60/212671, filed June 20, 2000, the disclosure of which is hereby incorporated herein 
by reference. 

[002] The present invention relates, in general, to the field of volumetric analysis, and 

more particularly to flow titration utilizing feedback-based flow ramp reversal. 

[003] Titrimetry is one of the few classical analytical methods still in wide use, for the 

determination of both major and minor components (the latter most notably for the 
measurement of water by Karl Fisher titrations). Titrations are not limited to solutions. 
Indeed, the origin of titrimetry has been traced back to Geoffroy in 1729; he evaluated 
the quality of vinegar by noting the quantity of solid K 2 C0 3 that could be added before 
effervescence ceased. 

[004] The foundation of volumetric analysis, as it is presently known, was laid by Gay- 

Lussac between 1 824 and 1 832. Mohr is especially credited for popularizing volumetric 
analysis, through the 1855 publication of his classic treatise on titrimetry. Compared to 
competing techniques, titrimetry exhibits excellent precision, convenience and 
affordability, but it is generally confined to a batch operation with slow throughput and 
requires significant amounts of sample and titrant. As long as the titrant concentration 
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is exactly known and volumetric ware needs no further calibration, true titrations require 
no calibration curve. This can be important in situations where it is difficult to prepare 
a pure standard solution of the analyte being titrated. 

[005] In routine analytical laboratories, automated titrators, often coupled to robotic 

workstations, have largely replaced manual titration. The hardware and algorithms for 
end point detection and the consequent feedback to control the rate of titrant addition 
have become increasingly sophisticated. In the research laboratory, titrators have been 
greatly miniaturized and automated and titrations have been demonstrated down to the 
femtoliter scale. The intrinsic batchwise nature of this measurement has however, 
remained unaltered in such efforts. 

[006] More than three decades ago, the paradigm of titrations from the volume to the 

flow domain was introduced. The basic arrangement involves two pumps; one pumps 
the sample at a constant rate while the other, a feedback based servo-controlled pump, 
delivers the titrant. The two streams are merged, and some desired property of the 
mixed stream (such as electrode potential) is measured by a suitable detector. In 
response to the detector output, the titrant pump changes speed and attempts to maintain 
the detector output at some constant level. This general system has been copied, 
reinvented, reconfigured, and (over)simplified numerous times since its original 
description, and linear titration plots have been shown to be possible. 
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However, manual rather than feedback based control of the titrant pump is still 
being advocated, with finite pauses at discrete mixing ratios so as to allow the detector 
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output to reach the equilibrium value reflective of that mixing ratio. This generates a 
complete multipoint titration curve. Such seemingly retrograde movements highlight 
the fact that the same problem faced by batch mode autotitrators reappears in the flow 
domain, perhaps in a more aggravated fashion. In the batch titration mode, following 
incremental titrant addition, a finite mixing and detection time is necessary before a 
detector reading is meaningful. In flow titrations, aside from the above factors, the time 
for the fluid transport to the detector is additionally necessary after a flow rate change. 
Altogether, these contribute to the overall lag time between the equivalence 
composition (or some other desired set point) being first reached at the sample-titrant 
confluence point and this actually being registered by the detector. Mixing is a more 
problematic issue in flow systems because it is difficult to incorporate active mixing 
(magnetic stirrers etc.) without incurring additional mixing volume that increases the lag 
time further. Detector response time can also be factor. Optical detector response is 
rapid but potentiometric detection, especially with a glass electrode, can be slow. 
Because of the sigmoidal nature of most titration curves (near the end point the change 
in the detector output is most rapid), the existence of a finite lag time is particularly 
important. Autotitrators solve this problem by reducing titrant delivery speed in 
proportion to the rate of change of the detector output, even performing discrete small 
incremental additions near the end point for highest accuracy. There is no barrier to 
implementing the same methodology in the flow domain, to perform measurements with 
discrete changes in the flow ratio while waiting at each step for a steady state reading 
and with progressively smaller changes as the end point is approached. Either this 
approach or the concept of holding the system at the end point in a steady state under 



feedback control has the virtue that they are true titrations. As long as the pumps are 
calibrated (in much the same way that volumetric ware used in conventional titrations 
have implicit calibrations), no further system calibration with samples of known 
concentration is needed. They have, however, the considerable disadvantages of a 
substantial time (>5 min) and sample/titrant consumption per titration. For applications 
as on-line process titrators, the significant time involved can cause a real problem in 
streams with fast-changing compositions. 

[008] Many ingenious and innovative approaches to solving these problems have 

emerged. The detector output has been used as the index to perform a half-interval 
search in the sample:titrant ratio to attain the equivalence point with excellent precision 
in under 3 min per titration. With air bubbles at each end of a trial mixture to isolate it 
from the next mixture (monosegmented flow), the method could be used for even slow 
detectors without a major increase in the time spent. 

[009] In such a binary approach to attain the end point, a maximum of n compositional 

trials are necessary for a desired accuracy of 1 part in 2 n . One approach is to use two 
variable speed pumps, which alternate their roles as sample/titrant pumps with the help 
of two 3-way valves. A titration is performed one way and then the pumps change their 
role to perform the titration again and an average of the two compensates for the lag 
time. Further, the temporal flow profile for the pumps are exponential, rather than 
linear, to provide better accuracy. 



[001 0] There are also other techniques, even though they are not true titrations in that 
calibration is required, such as concentration gradients in the titrant, rather than a flow 
gradient; such gradients can be generated by exponential mixing or by electrochemical 
generation of the titrant. 

[0011] In more recent years, flow injection titrations have been introduced. In this 
approach, the width of an injected response peak at some set point, often the 
equivalence point, is measured and is linearly proportional to the logarithm of the 
injected concentration. Although gradient dilution chambers are advocated, they are 
not essential. While logarithmic response permits a large dynamic range, precision is 
limited. To get 1 % accuracy, analysis time can be several min per sample and changes 
in sample viscosity can compromise accuracy. Instead, using the configuration of 
Figure 1 , the time to reach a titration end point, beginning from a zero titrant flow, as 
deduced from an indicator color change occurring as a result of a linear change in 
titrant flow rate, is linearly related to the sample concentration. This may be a more 
attractive approach in terms of time and precision. Further, for the purist, it is more 
likely to qualify as a titration. However, the problems of reducing the time for carrying 
out an accurate titration still remains. 

SUMMARY OF THE INVENTION 
[001 2] The present invention is directed to a new paradigm for continuous flow titrations 
by feedback-based flow ratiometry in which the lag time between sample-titrant 
confluence and detection of the status of the titration is made constant. The error in 
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measuring the titration results due to a lag time, whether the lag is due to the resident 
time in the mixer or is due to the detector response time, is continuously compensated 
for by averaging rapid backward and forward titrations. 

[0013] In accordance with the invention, continuous on-line titrations are based on 
feedback-controlled flow-ratiometry where the ratio of sample flow to titrant flow is held 
equal to the ratio of sample concentration to titrant concentration, and on the principle 
of compensating errors. System and methods operating under these principles have 
been thoroughly tested by applying them to acid-base neutralization titrations with 
indicator-based end-point detection, and will be described in terms of such tests for 
convenience. However, it will be understood that the invention is not limited to the 
titration of these particular samples or titrants. In a typical case, a total flow F T , 
consisting of sample and titrant flows, is held constant while the titrant flow F B varies 
linearly in response to a controller output voltage. The sample (e.g., an acidic solution 
to be titrated) flow F A constitutes the makeup (the difference between the titrant flow F B 
and the total flow F T ) and is added to the titrant at a point of confluence to provide the 
mixed stream F T Thus when the titrant flow is varied the sample flow also varies in the 
reverse direction (F A = F x - F B ). The status of the mixed stream F T is monitored by a 
detector and used either for governing the controller output or for interpreting the 
results of the titration. The titrant may be a standard base solution containing an 
indicator, in which case the detector may be an optical detector. However, it will be 
understood that other titrant solutions may be used to vary the properties of the mixed 
stream, with the detector being appropriate to the property to be detected. 
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[0014] Three methods of control were tested: a PID based control, a fixed flow rate 
control and a feedback-based flow rate control, with the last-named method being the 
preferred embodiment of this invention. In this latter approach, the titrant flow is initially 
ramped upward in accordance with a preselected flow rate pattern. At the instant a 
change in a selected property of the mixed stream is sensed by the detector, the actual 
titrant flow rate F H (which is produced by the upwardly ramping flow control signal) is 
higher than the true equivalence flow rate F E because of the lag time between the 
occurrence of the first property change and its detection. The sensing of the change 
in property is used to cause the system controller output to immediately reverse its 
ramp direction so that the titrant flow is ramped downwardly in accordance with the 
same flow rate pattern. At the instant another change in property, in the opposite 
direction this time, is detected, the titrant flow rate F L (produced by the downwardly 
ramping flow rate control signal) is lower than F E by exactly the same amount that F H 
was higher than F E . This is the principle of compensating errors (F E = (F H + F L )/2) 
which allows true titrations with excellent reproducibility and speed (0.6% rsd @ 3 
s/titration and 0.2% rsd @10 s/titration), and reduced titrant volume consumption (as 
little as 12 uL/titration), and which solves an old conceptual problem in flow-based 
titration. 

BRIEF DESCRIPTION OF DRAWINGS 

[0015] The foregoing, and additional objects, features and advantages of the present 
invention will become apparent to those of skill in the art from the following detailed 
description of preferred embodiments, taken with the accompanying drawings, in which: 
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[0016] Fig. 1 is a diagrammatic illustration of a number of prior titration systems; 

[0017] Fig. 2 is a diagrammatic illustration of a titration system in accordance with the 

present invention; 

[0018] Fig. 3 is a diagrammatic illustration of a system for varying the sample 

concentration in the system of Fig. 2; 
[0019] Fig. 4 illustrates controller output vs time for the system of Fig 3 using a PI D flow 

rate controller; 

[0020] Fig 5. Illustrates controller output vs time for a triangular wave flow rate 
controller, showing detector output D out ; 

[0021] Fig. 6 illustrates, in graphs (a) and (b) the detector output vs controller voltage 
for a triangular wave flow rate controller over 10 titrations each of a sample of (a) 50 
mM and (b) 100 mM of HCI; 

[0022] Fig. 7 illustrates, in graphical form, the voltage output vs time for a feedback- 

based controller in accordance with the present invention, having a triangular flow rate 
control voltage; 

[0023] Fig. 8 is a graphical illustration of the feedback-based triangular flow rate control 
method of the invention utilizing a continuous change in analyte concentration, 
illustrating in graph (a) function generator output FG 0Ut that governs the analyte 
concentration and controller output V c vs. time, and illustrating in graph (b) the 
reciprocal of the controller output at equivalence (1/V E ), which tracks the function 
generator output; and 

[0024] Fig. 9 is a graphical illustration of the controller output vs detector output for two 
different analyte concentrations in a feedback-based triangular wave control method 



in accordance with one embodiment of the invention. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
[0025] As previously described, Fig. 1 illustrates the known Blaedel-Laessig titration 
configuration 10, wherein a sample S is supplied by a pump 12 to a mixing coil MC 
through a supply line 14. The titrant T is supplied to mixing coil MC through a variable 
pump 16 and line 18, and a detector D measures property changes in a mixed stream 
in flow line 20 from the mixing coil, with the waste fluid W flowing from the detector by 
way of outlet line 22. 

[0026] The basic flow diagram of a titration system 30 in accordance with the present 
invention is shown in Fig. 2, wherein a variable titrant pump 32 supplies a titrant T 
through an output line 34 to produce a flow F B to a mixing reactor MR. A sample 
source S supplies a sample through line 36 to produce flow F A , with the titrant flow on 
line 34 aspirating the sample at junction 38 and carrying it to the mixing reactor MR. 
In prior configurations, variations in pump flow rates resulted in variations in total flow 
rate and a fixed hardware arrangement produced corresponding variations in the lag 
time between sample titrant confluence at function 38 and subsequent detection. To 
overcome this variation, the configuration of Fig. 2 includes a pump 40 which keeps the 
total flow rate F T from the reactor MR constant, and includes a flow rate controller (to 
be described) for titrant pump 32. The flow rate F B of the pump 32 is adjusted to be 
equal to or just below F T . The sample flow rate F A represents the difference between 



F T and the titrant flow rate F B . Since the mixed stream F T from pump 40 to reactor MR 
on line 42 varies between 100% sample and 100% titrant, it is possible in principle to 
titrate a sample of any concentration with a titrant of any concentration. However, to 
obtain good precision and accuracy, a judicious choice of titrant concentration is 
appropriate, based on the sample concentration. 

[0027] In experiments, the sample concentration was varied as a function of time using 
the arrangement illustrated in Fig. 3. A ramp generator 50, such as a Tektronix FG 
504 function generator, was connected to a voltage-controlled pump 52 to vary the flow 
of water in line 54 to dilute a constant flow of a sample stream on line 56 provided by 
a constant rate pump 57. The water and the sample were mixed at a mixing coil MC 
to provide a mixed stream on line 58. Part of this mixed stream on line 58 was 
aspirated through the sample aspiration line 36 shown in Figure 2, while the rest was 
allowed to go to waste (W). 

[0028] Variable speed pumps such as Gilson Minipuls 2 or Rainin Rabbit- 
Plus/Dynamax pumps, having 10 stainless steel rollers, were used for the pumping 
needs. The variable pumps 16, 32 and 52 may be externally voltage-controlled with a 
0-5 V DC analog input. For good mixing without significant residence time, single bead 
string reactors MR^ and MR 2 , each incorporating, for example, a single strand chain of 
beads with an average bead diameter of 0.5 mm, were contained in a tube having 0.81 
i.d. and 25 mm long, were used on both the inlet and outlet of the final pump 40, 
illustrated in Fig. 2. Residence time for mixing is less critical in the system of Fig. 3, 
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so a mixing coil MC (600 mm x 0.66 mm i.d.) indicated at 59, suffices for this purpose. 
Except as stated, 0.51 mm i.d. Pharmed pump tubes, or lines, were used throughout, 
with the total flow rate F T held constant at 1 .9 mL/min. 

| Although changes in various properties of the total flow may be measured to 
determine titration end points, the use of indicators in the titrant to change the color of 
the sample is common, and the following description of tests of the present invention 
will be in terms of such indicators in combination with optical detectors. It will be 
understood, however, that the invention is not limited to such indicators or to optical 
detection methods, but is equally applicable to, for example, potentiometric or 
conductometric detection methods. In one example of the invention, a conventional 
indicator was premixed in the titrant T and the optical absorption of the indicator 
present in the mixed stream F T was measured by a simple on-tube Light Emitting Diode 
(LED)-photodiode (PD) based detector 62. The detector 62 used in the experiment 
described above consisted of a 1/4-28 threaded male-male union with a center partition 
made for chromatography (P/N 39056, Dionex) with an LED 64 emitting at 605 nm (P/N 
HAA5566X, Stanley Electric, Tokyo) on one side and a silicon photodiode 66 (PD, P/N 
BPW 34, Siemens) on the other, transparent tubing 68 of FEP Teflon (0.8 mm i.d., 1 .2 
mm o.d.) carrying the mixed stream F T passed between LED 64 and PD 66 in a 
perpendicular fashion. The PD produced an output 70 which was supplied to a current 
amplifier 72 (Amp, Model 427, Keithley) which typically was set at its minimum 
response time of 1 0ps, and the resulting amplifier output voltage on line 74, which was 



linearly related to the optical transmittance of tubing 68 and its contents, was supplied 
to a controller 76. 

[0030] Two different controllers were utilized in the above-described experiments. The 
first controller was a commercially available PID type process controller (Omega 
CN761 60), while the second incorporated a personal computer (PC)-based system with 
a control algorithm written as described below. As an intermediate step in operating the 
latter system, a function generator (Tektronix FG 504) was used in place of the PC 
based system to ramp the flow produced by the titrant pump 32 up and down, in a blind 
fashion. The output from the amplifier 72 and from controller 76, as well as other 
operating parameters were acquired on a PC using a 12-bit data acquisition card 
(DAS-1601, Keith ley/Metrabyte). In the case of the fully PC-based system, a PCMCIA 
card (PCM-DAS16D/12AO 12-bit A/D, D/A, Computerboards Inc., Middleboro, MA) 
housed in a Pentium II- class notebook computer (Latitude, Dell Computer Corp.) was 
used. Commercially available reagents of analytical reagent grade (indicator grade for 
indicators) were used without further purification in the task of the present invention. 

[0031] Tests were carried out using the configuration of Fig. 2. In this configuration, at 
the equivalence point of the titration the following equation will hold: 

C A (F T -F E ) = C B F E (Eq. 1) 

where C A and C B represent the concentration of sample and titrant respectively, and F T 
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is the invariant total flow rate of the system. The titrant flow rate F B is designated F E at 
the equivalence point. In these tests, generally the sample was an acid and the titrant 
was a standard base solution already containing an indicator. The value of F B was 
linearly related to the output voltage V c from the controller 76. Thus, Eq. 1 is rewritten 
as follows: 

C A (F T -KV E ) = C B kV E (Eq.2) 

where k is a constant of proportionality, kV E is equal to F E , and V E is the value of V c at 
the equivalence point. Therefore, 1/V E is proportional to 1/C A in the configuration 
shown in Fig. 2, all other terms in the following equation being constant: 

(V E )- 1 = (kC B /F T )(C A )" 1 + k/F T (Eq. 3) 

[0032] Proportional Integral Derivative (PID) controllers are widely used for control of 

temperature, pressure and other process parameters, and many controllers of this type 
are commercially available. In a test, such a controller was connected to maintain the 
system of Fig. 2 at such a titrant flow rate that the mixed stream F T was exactly 
neutralized in a strong acid-strong base titration. The titrant was 100 mM NaOH 
containing 0.2 mM bromthymol blue (BTB); sample: 50-200 mM HCI. There was no 
buffer capacity at the equivalence point except for that provided by the indicator, and 
a very slight deviation caused an indicator color change. The detector output at 74 was 
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measured with the mixed stream distinctly acidic, where the indicator was completely 
yellow, and distinctly basic where the indicator was completely blue, and the controller 
was set to maintain the detector output at the midpoint of these acidic and basic values. 
Under these conditions, if steady state control can be achieved, V c will directly reflect 
analyte concentration since F B is linearly related to V c . 

[0033] Self-tuning abilities are provided in most modern microprocessor based PID 
controllers to find optimum values of P (gain), I (bias) and D (time constant) to maintain 
good control. For the present system, where the detector output is practically bistable 
with a sharp transition between two states (indicator blue and indicator yellow), self- 
tuning was ineffective to maintain control. The controller output oscillated between its 
highest and lowest permitted values. Manual settings of the PID values led to better 
results; the control voltages for two runs is shown in Fig. 4 at graphs 80 and 82. In 
graph 80, the PID values were set at P = 25, I = 0.1 and D = off, and the controller 
output was obtained as the HCI concentration in the sample was varied from 50 mM to 
200 mM and back to 5 mM in seven steps, illustrated at points 1a to 1d. In the test, 
step 1a = 50 mM, step 1b = 100 mM, step 1c = 150 mM, and step 1d = 200 mM. In 
graph 82, the same steps in sample concentration were used, at 2a - 2d, with the PID 
values set at P = 100, I = 0.1 and D = 0.01. 

[0034] It was found that some degree of control was possible with high P and low D 

values. With an increase of P value, oscillation became less significant, but it took a 
longer time for the system to stabilize at the set point. No significant improvements 
were observed by using greater time constants for the detector to reduce detector 
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noise, or by reducing lag time by placing the detector on the aspiration side of the final 
pump 40 (Fig. 2). It was also found that the degree of oscillation for the same sample 
concentration can be unpredictable, as illustrated by graphs 80 and 82. More 
importantly, if the midpoint of each oscillation is V E , the relative standard deviations 
would be significant. Accordingly, PID control is not particularly well suited for 
maintaining a system at the equivalence point where the rate of change is very steep 
near the set point. 

[0035] The foregoing difficulty can be overcome, in accordance with the preferred form 
of the present invention, by scanning the titrant flow in the vicinity of the equivalence 
point, without attempting to keep the mixed effluent at equivalence. If the titrant flow is 
being ramped upwards linearly, for example, at the instant a change in the color is 
sensed by the detector, the titrant flow rate F H is higher than the true equivalence flow 
rate F E because of the lag between the time equivalence is reached in the mixing reactor 
MR and the time the mixed liquid reaches the detector, for during that time lag, the titrant 
flow continues to increase. Calibration of an entire system can include implicitly taking 
into account this lag time. F E can also be expressed as: 

F E = kV E = F H - rt lag (Eq. 4) 

Where r is the ramp rate (dF B /dt) and t, ag , the lag time, is a combination of the physical 
transit time from the confluence point and the detector response time. There may be 
nothing wrong with calibrating the system with standards for subsequent assays except 
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that it ceases to be a true titration; calibrations beyond flow rates are required. This may 
be a nontrivial issue; if a system is calibrated for use with one set point, it may have to 
be calibrated again for use with another set point, since for example, the precise 
response times of a glass electrode varies depending on the pH regime. 
[0036] If, however, the titrant flow F B is then decreased from its high value at the same 

ramp rate r, and the instant of a change in color, in the opposite direction this time, is 
again sensed, the measured titrant flow rate F L will be lower than F E , in a mirror image 
fashion of the previous situation: 

F E = kV E = F L + rt lag (Eq. 5) 

F E may then be calculated from Eqs. 4 and 5 without the use of system calibrations and 
without knowing the specific values of r and t lag , thus permitting true titrations: 

F E = <F H + F L )/2 (Eq. 6) 



[0037] We refer to this herein as the principle of compensating errors. 

[0038] This concept was tested using a controller 76 that is simply a generator producing 
a triangular wave V c , where the triangular wave is used to control F B , as illustrated at 
curve 84 in Fig. 5. Actual data output, utilizing the sample titrant described above, with 
a 100 mM HCI, at a ramp rate of 100 mV/s is shown in Fig. 5. In response to the 
triangular wave V c controlling the titrant flow, the detector output D out on line 74, 
illustrated at 86 in Fig. 5, basically executes a rectangular wave pattern. The yellow 
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form of the indicator had practically no absorption at the monitoring wavelength. Thus 
at the high end, D out is flat. On the other hand, even after the indicator turns blue, further 
increase of F B brings still more indicator in the system, since the indicator is 
incorporated in the titrant, and at the low end, D out executes a shallow V, with the bottom 
of the V being approximately temporally coincident with the apex of V c , the difference 
being the lag time, t, ag . 

[0039] The effect of the ramp rate r was examined at constant scan limits (1 .0 to 4.5 V, 
20% - 90% of maximum pump rate) using the same NaOH-HCI system as in the previous 
paragraph. Pump flow rates were calibrated so that the correct value of F E (or V E ) would 
be known a priori. At low to moderate values of r (0.0302, 0.0607 and 0.1203 V/s), V E 
values obtained were in excellent agreement with the expected true value and exhibited 
good precision, with a relative standard deviation (RSD ) of less than 1%). However, 
at high scan rates (0.24 V/s), the experimental V E values were inaccurate. Under these 
conditions, and with the t, ag in the system used for these tests, V c reached its limits 
before a transition in D out was detected. 

[0040] Since the number of titrations that can be conducted per unit time increases with 

increasing the ramp rate, or scan rate, r, further experiments were conducted with r = 
0.100 V/s. HCI solutions (50, 75, 125, 150 and 200 mM) as the sample were measured 
with fixed scan limits of 1.0 to 5.0 V. Plotting the resulting data (1/V E vs. 1/CHCI) 
exhibited good linearity (linear r 2 0.9976). A plot of D out as a function of V c is shown in 
Fig. 6 for 10 titrations each of (a) 50 mM (graph 88) and (b) 1 00 mM HCI (graph 90). V E 
essentially represents the abscissa value corresponding to the center of mass of the 
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parallelogram of each graph. Compared to extant literature methods, this approach is 
quite competitive (80 s/cycle, <1% precision), but Figure 6 also shows very clearly that 
the scheme results in large amounts of time being spent in a useless manner. For 
example, in graph 88, the system unnecessarily scans in the VC 3 to 5 V range and 
similarly it spends unnecessary time in the titration of graph 90 in the VC, to 3 V range. 
[0041] A more efficient method than the one described above, and the preferred 
embodiment of the invention, involves reversing the direction of the titrant pump as soon 
as the equivalence point is crossed. This is accomplished by sensing the detector 
output and changing the direction of the control voltage V c as soon as some preset 
threshold in D out is crossed. Other standard endpoint functions based on the first or 
second derivatives of D out can also be implemented, if desired. Both the principle and 
the results are illustrated in Fig. 7, where graph 92 is the controller 76 output V c and 
graph 94 is the detector voltage D out to the controller 76, by using the same titrant and 
sample as in Fig. 5. In contrast to the previous method, controller scan limits are not 
fixed but are, in effect, V H and V L and thus vary with the concentration of the analyte. 
As soon as D out indicates an alkaline mixture, where V H has been reached, the controller 
output V c is ramped downward. When D out indicates an acidic mixture, where V L has 
been reached, controller output V c is immediately ramped upward. The resulting V c 
waveform 92 has a constant frequency of (4t, ag )" 1 and displays an amplitude of 2rt !ag . 
Note that these properties of the V c waveform are independent of the analyte 
concentration. 



[0042] The frequency and amplitude properties of the controller output waveform 
suggest the possibilities of diagnosing and/or compensating for flow inconstancies, 
whereas the DC component of the controller output V c is related to the analyte 
concentration. This DC bias moves up or down as the concentration of the analyte 
increases or decreases. The detector setpoint (D out , high) at which V c begins a 
downward ramp does not have to be the same as the detector setpoint (D out , low) at 
which the ramp goes back up. This is of practical importance since all real signals 
contain some noise. When these points are set identically, false triggering, such as 
premature ramp reversal in either direction, can and will occur. To avoid such problems, 
D out , high, should differ from D out , low, by at least 2 times the detector noise. Because 
the transition is very steep, it makes no real difference in the ultimate results in V E . 

[0043] For a discrete sample, at least one V H and one V L value are necessary to 

compute V E . However, when the system is being applied continuously to a flowing 
stream, V E values will be computed by averaging the most recent V H or V L value with the 
immediately preceding V L or V H value. Since the period of the V c waveform is directly 
dependent on the lag time of the system, t| ag , it is essential to reduce it to improve 
throughput, but t iag cannot be reduced indefinitely without affecting the completeness of 
mixing and thus increasing detector noise and decreasing system reliability. These 
interrelated issues are of critical importance. 

[0044] The effects of positioning the detector 62 upstream or downstream of the final 
pump 40 and of the type and the size of the reactors MF^ and MR 2 (for example, knotted 
tubing, single bead string reactor, or a porous frit-tee) at various scan rates and sample 



(HCI)concentrations were examined, and a representative set of results is presented inTable 1. 



Table 1. 


Effect of detector position and reactor type 




Detector Position 


D a 


D 


D 


D 


D 


U b 


U 


Reactor 


none 


KT C 


KT 


KT 


two SBSR d 


SBSR 


SBSR 


Reactor length, cm 


NA e 


15.0 


30.0 


60.0 


2 x2.5 


2.5 


10.0 


Titration time, s 


10.4 


12.6 


15.5 


20.1 


10.2 


3.2 


4.6 


% RSD of V E , (n=20) 


0.32 


0.46 


0.40 


0.35 


0.22 


0.56 


0.76 



: a Downstream of final pump. b Upstream of final pump. 0 Knotted tubing, 0.81 mm i d . 
d Single bead string reactor. e Not Applicable. 
sfTitrant: 100 mM NaOH containing mM BTB, sample 100 mM HCI, r = 0.100V/S 

[0045] It was possible to reduce titration time to as little as 3.2 s with a small penalty to 
the precision in V E (RSD 0.56 %) by locating the detector upstream of pump 40 (Fig. 
3). This loss of precision (although acceptable for many purposes) was attributable to 
increased detector noise since the active mixing provided by the pump tubing was no 
longer available. To improve pump induced mixing, a combination of high pump 
rotation rates and small bore pump tubing, rather than large bore tubing and lower 
rotation rates, was chosen for the desired flow rate regime. A porous frit (pore size 90- 
1 30pm) - tee proved to be incompatible with peristaltic pumps because of its substantial 
resistance to flow. 

[0046] The effect of changing the ramp rate r for controller 76 was examined over the 
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0.010 - 0.200 V/s range with the same titrant-sample combination as used 
Table II provides a summary of the results. 



in Table I. 



Table II. Effect of Scan Rate 



% ; Titration 



an rate, V/s 


I v H 


L.SDvH a _j 


v L 


SD VL a 


v E 


j SD VE a 


I RSD VE 


Time, s 


0.01 


i 3.49 


i 0.01 i 


3.35 


i 0.02 } 


3.42 


i 0.01 


! 0.40 


14.3 


0.025 


i 3.52 


0 02 


3.24 


i 0.02 ' 


! 3.38 


i 0.01 


I 0.43 


11 6 


0.05 


i 3.68 


0 02 


3.14 


0 01 


: 3.41 


i 0.01 


I 0.35 


108 


0.075 


I 3.80 


i 0.01 i 


3.03 


i o.oi ; 


3.41 


i 0.01 


i 0.26 


10 35 


0.1 


I 3.94 


0 01 


2.92 


! 0.01 i 


; 3.43 


i 0.01 


i 0.23 


i 10.2 


0.15 


I 4.16 


0 03 


2.70 


i 0.02 ; 


i 3.43 


i 0.02 


i 0.48 


9.8 


0.2 


4 39 


0 02 


2.48 


! 0.03 I 


i 3.44 


! 0.02 


i 0.44 


9.6 



: a Standard deviations of preceding parameters. Wherever applicable, the units are volts. 



[0047] In principle, the scan rate (or ramp rate) r should not have a direct influence on 
the titration time. This was at least approximately true; a 20-fold increase in r resulted 
in only a 33% decrease of the titration time. The limited effect of r on the titration time 
that was observed is a practical consequence of a finite mixing and detection volume. 
The V E values were virtually constant irrespective of the scan rate and the observed 
range at different scan rates was within 0.7% of the mean, and this range included the 
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independently determined true value. V H and V L increasingly diverge from V E as r 
increases, in accordance with equations. 4 and 5. A plot of V H vs. r should thus have 
V E as the intercept and a slope equal to t, ag . The data in Table II yields a value of V E 
= 3.430 V, a lag time t lag = 4.87s, and a linear r 2 value of 0.9979. Similarly, a plot of V L 
vs. r exhibits a linear r 2 value of 0.9979, an extrapolated V E of 3.366 V, and yields t lag 
= 4.47s. This range of V E is almost the same as the range of V E values observed by 
numerical averaging of successive V H and V L values at individual scan rates, lending 
additional credibility to the theoretical basis of the present invention. 

[0048] Although the scan rate was not a major factor in determining the titration time for 
constant or very slowly changing sample concentrations, a higher scan rate will reach 
V H or V L values more rapidly when the sample concentration changes quickly. 
However, a very fast scan may result in V c reaching the scan limits before a transition 
in D out is detected (as was noted above for a function generator controlled operation of 
the system without feedback, for r > 0.24 V/s.) In addition, a slower scan may improve 
the accuracy and precision in determining V E . In looking at the % RSD values for V E , 
the RSD values increase at r > 0.100 V/s and also increase at r < 0.05 V/s. It was 
concluded that the 0.05-0.100 V/s, which corresponds to 1%-2%/s of the maximum 
possible pumping rate, is the best range in which to operate. 

[0049] Although the herein-described experiments utilize a linear scan rate r produced 
by a triangular wave output from controller 76, it will be understood that not only can 
this scan rate be changed, it can also be varied during the scan to produce control 
voltages having different waveforms, with different increasing and decreasing slopes. 
Thus, the invention is not limited to a particular controller output pattern, or wave shape. 
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[0050] For virtually any process application, the sample consumption rate is not 

particularly important, whereas minimizing the titrant consumption is highly desirable 
since this decreases the replacement frequency of a reagent that must be carefully 
made. The utility of more concentrated NaOH solutions as titrant, in conjunction with 
smaller diameter (0.25 and 0.44 mm) pump tubes in addition to the original 0.51 mm 
diameter tube used in all the above experiments was, therefore, examined. Typical 
results are presented in Table III. 

Table III. Effect of Titrant (NaOH) Concentrate 3 



Tubing i Nominal i i t i i % t Titration? 



i.d. 
mm 


! C B . M 


I v., 


SD VH ° 


v L 


SD VL 


v E 


SD VE 


RSD VE 


ii Time, s 


0.25 


i 0.5 


I 4.65 


j 0.0148 \ 


3.62 


0 0170 


4.14 


00112 


0.27 


\ 10.35 


0.25 


I 1.0 


1 2 76 


0 0115 


1.70 


i 0.0116 i 


2.23 


0 0081 


0.36 


1 10.55 


0.25 


| 2.5 


1 43 


) 0.0187 j 


0.37 


0 0134 


0.90 


00115 


1.28 


i 10.60 


0.44 


? 0.5 


2.02 


0 0095 


0.96 


0 0084 


1.49 


0 0063 


; 0.42 


| 10.50 


0.44 


| 1.0 


1 31 


| 0.0071 j 


0.27 


0 0091 


0.79 


j 0.0057 I 


0.72 


| 10.40 


0.51 


0.1 


3 81 


0 0084 


2.80 


0 0132 


3.30 


0 0078 


0.24 


| 10.10 


0.51 


0 


1 57 


0 0130 


0.58 


I 0.0112 i 


1.07 


0 0085 


0.79 


9 90 


0.51 


! i.o 


I 1.01 


0 0055 


0.05 


0 0098 


0.56 


0 0056 


1.01 


j 10.15 



i: a 0.1 M HCI is sample, r = 0.100 V/s. 

§n b Standard deviations of preceding parameters. Wherever applicable, the units are Volts. 

[0051] A combination of low C B and narrow tubing is insufficient to titrate the sample 
and was not tested. Likewise, a combination of high C B and wider bore tubes is 
impractical. In general, the absolute standard deviations decreased with an increase 
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in NaOH concentration. However, not surprisingly, the relative standard deviation 
increased with increasing C B because V E decreased more significantly compared to the 
standard deviation (SDV E ) in V E . Nevertheless, the use of a small diameter (0.25 mm 
i.d.) tube allowed the relative standard deviation of V E (RSDVE) for titrations using 500 
and 1000 mM NaOH to be kept well below 0.5%. Similar results, not shown, were 
obtained for a scan rate of 0.050 V/s. By using 1000 mM NaOH instead of 100 mM 
NaOH, the titrant consumption could be reduced to 18 % of the original value. 
[0052] Note that in the configuration of Fig. 2, the average volumetric consumption of 
the titrant (which is assumed to be linearly related to F E ) is linearly related to 1/(1 + 
C B /C A ). (See Equation 1). While large gains are made initially in F B consumption with 
increasing C B , there are diminishing returns on an absolute scale at higher and higher 
titrant concentrations. Nevertheless, it is remarkable that with a titrant concentration 
25 times that of the sample, it is still possible to perform titrations with a precision only 
slightly over 1%, at only -10 s/titration and consuming 1 1 .7 uL/titration. The precision 
at very high titrant/sample concentration ratios may be improved by using even smaller 
diameter pump tubes. 

[0053] To monitor compositional changes in a stream accurately using this technique, 
it must be possible to carry out each forward and backward titration at a much faster 
rate than the rate at which the stream is changing composition. 

[0054] A process stream in which the reciprocal of analyte (HCI) concentration changed 
linearly with time was created by the system depicted in Fig. 3. The dilution flow to a 
constant flow stream of an acidic analyte was increased linearly with time by a slow 
triangular wave output (FG out ) of function generator 50. With a FG out cycle time of 93.67 



min, the minimum and maximum analyte concentrations were 50 and 180 mM. Since 
the entire range is spanned within one half cycle, this means the analyte concentration 
varied by almost a factor of four in a period of ~45 min. This degree of change more 
than adequately represents the maximum change that occurs for a critical and major 
component in a real process stream. Six hundred and eight titrations were made within 
one cycle of the function generator 50, resulting in an average titration time 
requirement of <9.2 s per measurement. Graph (a) of Fig. 8 shows V c (curve 94) and 
FG out as a function of time. The reciprocal of V E is shown similarly as a function of time 
in graph (b) of Fig. 8 at curve 98, and this linearly tracks FG out as would be expected 
from theory. 

[0055] The system of the present invention was applied to a number of acid-base 
neutralization titrations in addition to HCI-NaOH: for example CH 3 COOH-NaOH, H 3 P0 4 - 
NaOH and NH 3 (aq)-HCI. Indicators were selected that not only have a pK, n (indicator 
dissociation constant) value in the desired range but also that are blue in one form 
(basic) so that the 605 nm LED detector 62 could be used without any further 
modification. For H 3 P0 4 , titrations both at the first and second equivalence points 
were carried out using a separate indicator for each. The results are summarized in 
Table IV. 
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Table IV. Results for Different Acid-Base Neutralization Reactions 



Number of 
Individual 



i Sample 


Titrant 3 


i Indicator 


i Scan Rate 1 
* V/s 


i Measurement ? 
range, mM 


Concerns 
Examined 


1 Linear 

3 r2 


nOI 


NaOH 


BTB b 


3 0.05 


50 - 200 


6 


1 0.9989 


HCI 


NaOH 


BTB 


0 1 


50 - 200 


6 


I 0.9991 


CH3COOH 


NaOH 


TB° 


i 0.05 


25 - 200 


7 


0.9994 


CH3COOH 


NaOH 


TB 


0 1 


25 - 200 


7 


0 9994 


H 3 P0 4 


NaOH 


BCG d 


0.05 


12 5 -200 


9 


i 1.0000 


H 3 P0 4 


NaOH 


BCG 


0 1 


! 12.5-200 ! 


9 


I 0.9996 


H 3 P0 4 


NaOH 


jpe 


I 0.05 


12 5- 125 I 


9 


1 0000 


H 3 P0 4 


NaOH 


i TP 


1 0.1 


12 5- 100 I 


9 


! 1.0000 


NH 3 (aq) 


HCI 


BCG 


0 05 


25 - 200 


7 


1 0.9999 




HCI 


I BCG 


0 1 


25 - 200 


7 


I 0.9999 



ii a 100 mM contg. 200 uM indicator. b Bromothymol Blue. 0 Thymol Blue. 
;ii d Bromocresol Green. e Thymolphthalein 



[0056] The linearity refers to a plot of 1A/ E vs. the reciprocal of the analyte 
concentration. 

[0057] Illustrative V c vs. D out plots 100 and 102 for two different HCI 
concentrations are shown in Fig 9, with the data set for each plot representing a 
total of 20 titrations. A comparison with Figure 6 clearly indicates the substantial 



superiority of the feedback-based approach. Note that these titrations were 
conducted with identical upward and downward ramp reversal set points. 
Considering that these plots are very data dense (each plot contains-8500 actual 
plotted points and the data file size for each plot exceeds 1 MB), the number of 
errant points are remarkably few. 
[0058] The error compensated feedback based flow ratiometric titration method 
described here displays good precision coupled to unprecedented speed. It will 
be ideally coupled to continuously flowing streams, whether for measurement or 
control, in a variety of situations. It will permit the use of universal indicators and 
multiple wavelength detection for determining multiple analytes per titration for 
indicator based detection, and is applicable to other detection methods, such as 
detectors using pH electrodes. A pH electrode responds more slowly, increasing 
t, ag , but on an absolute scale, the measurement rate is still quite fast, requiring < 
15 s per titration. 

[0059] Although the present invention has been described in terms of preferred 
embodiments, various modifications and variations will be apparent to those of 
skill in the art. Accordingly, the true scope and spirit of the invention is limited 
only by the following claims. 
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